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Abstract 

The time-integrated CP asymmetry in the decay —)• K^Kg is measured using 
3fb“^ of proton-proton collision data collected by the LHCb experiment at centre- 
of-mass energies of 7 and 8TeV. The flavour of the meson is determined by use 
of the decay D*^ —?■ D^7r~^ and its charge conjugate mode. The result is 

Acp = -0.029 ± 0.052 ± 0.022 , 

where the first uncertainty is statistical and the second systematic. The result is 
consistent with Standard Model expectations and improves the uncertainty with 
respect to the only previous measurement of this quantity by more than a factor of 
three. 
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1 Introduction 


In the Standard Model, CP violation in charm decays is expected to be small and hence 
potentially sensitive to contribntions from New Physics. Although measurements of the 
time-integrated CP asymmetry in decays to pairs of charged mesons showed hints of 
CP asymmetry at the level of 0.7%, the combined results are not yet conclusive Of 

particular interest, both for the search of New Physics and for the understanding of penguin 
contributions, are decays of mesons into a pair of neutral mesons, such as the decay 

^ Kg Kg If the CP asymmetry in decays to charged mesons is conhrmed 

and assuming moderate breaking of the SU(3) flavour symmetry, the CP asymmetry of 
this mode could be of 0(1%) or even larger [^. From a more recent Standard Model 
based analysis of the contributing amplitudes, a 95% confidence level upper limit of 1.1 % 
for direct CP violation in the decay ^ Kg Kg has been derived j^. The single previous 
measurement gave Acp = (23 ± 19) % [^. 

Here, we present the hrst result from the LHCb collaboration on CP violation in the 
decay —)■ KgKg. The measurement is based on mesons originating from D*~^ —> D^7i~^ 
decays, where the flavour of the meson can be inferred from the charge of the “slow” 
pion from the D*~^ decay. Throughout this document, D*^ stands for Zl*(2010)“*’. As 
a control channel, the decay )■ K~7i^ is used to estimate production and detection 
asymmetries. The analysis uses 3fb~^ of proton-proton collision data collected with 
the LHCb detector in 2011 and in 2012, at centre-of-mass energies of 7TeV and 8TeV, 
respectively. 


2 Detector and simulation 


The LHCb detector [9,10 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing 5 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200GeV|3 The minimum distance of 
a track to a primary vertex (PV), the impact parameter, is measured with a resolution of 
(15 -|- 29/pt) M-ni, where px is the component of the momentum transverse to the beam, 
in GeV. The polarity of the dipole magnet is reversed regularly throughout the data-taking 
period, which allows to determine and correct for charge asymmetries due to the detector 
geometry. 

Different types of charged hadrons are distinguished using information from two 
ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identihed by a 


^The inclusion of charge-conjugate processes is implied throughout this article. 
^We use natural units where c = 1. 
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Table 1: Definition of candidate categories. 


Categories Description 

LL both Kg are of category long and not selected by the dedicated trigger 

LD one Kg is long, the other one is downstream 

DD both Kg are downstream 

LLtrig both Kg are of category long and selected by the dedicated trigger 


calorimeter system consisting of scintillating-pad and preshower detectors, an electromag¬ 
netic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of 
alternating layers of iron and multiwire proportional chambers. The online event selection 
is performed by a trigger, which consists of a hardware stage, based on information from 
the calorimeter and muon systems, followed by a software stage, which applies a full event 
reconstruction. 

In the simulation, pp collisions are generated using Pythia 11 with a specific LHCb 
configuration [T^. Decays of hadronic particles are described by EvtGen 13 , in which 


hnal state radiation is generated using Photos [^. The interaction of the generated 
particles with the detector, and its response, are implemented using the Geant4 toolkit 15 
as described in Ref. 


16 


3 Selection 

Signal decays are reconstructed in the decay mode D*~^ —with —)■ KgKg and 
Kg —)■ 7r'''7r“ [^. To collect as many —)■ KgKg decays as possible, events from all 
available physics triggers are considered. Candidate events are accepted if the four pions 
assigned to the Kg decays are sufficient to trigger the event or if the rest of the event, 
without the slow pion from the D* decay, satisfies a trigger condition. Excluding the slow 
pion from the trigger decision minimises any bias on the CP asymmetry due to the trigger. 

The Kg —)• 7r’''7r“ decays are reconstructed in two different categories: the first involves 
K! mesons that decay early enough for the daughter pions to be reconstructed in the 
vertex detector; the second contains Kg mesons that decay later such that daughter 
track segments are only reconstructed in the tracking detectors downstream of the vertex 
detector. These categories are referred to as long (L) and downstream (D), respectively. 
The less abundant long category has better momentum and vertex resolution than the 
downstream category. As the final state contains two Kg mesons, there are three possible 
combinations labeled LL, LD, and DD. A dedicated software trigger selection for the 
—)■ KgKg decay was implemented in 2012. As this trigger only accepts signal candidates 
composed of two long Kg, a fourth category, LLtrig, is dehned where the dedicated 
trigger accepted the signal candidate. The four categories are listed in Table [Tj 

The decay vertex of the candidate is reconstructed from the pion trajectories. 
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constraining the mass to its known valne 18 and the flight direction to point to 
a PV 17 . To rednce the contamination from the decay )■ iPgTr+vr”, only candidates 
with signihcant decay times for both decays are accepted. This requirement has a 
signal efficiency of over 99 %. The candidate is combined with a pion to produce a 
D* candidate. Fiducial cuts are applied on the kinematic properties of this slow pion to 
remove regions where the detection charge asymmetry is large. A cut on the invariant mass 
of the KgKg system of ±20 MeV around the known value of the mass is applied. The 
efficiency of this cut is 98% in all categories, except for DD, where it is 94%. Candidates 
are further selected by requiring the difference between the D*^ and candidate masses, 
Am = rriD* —moo, to be less than 155 MeV. 

To further reduce combinatorial background a multivariate analysis (MVA) method 
19,^ is used. It is based on a rule-based learner applying the methods of bagging and 
instance weighting (see, e.g., Refs. 


22,^). Separate samples have been used for training 


and testing, where simulated events have been used as signal proxy, while the background 
sample was a mixture of data from sidebands (the mass ranges 1764.84 — 1844.84 MeV 
and 1884.84 — 1964.84 MeV) and simulation. Kinematic quantities, decay time variables, 
geometric quantities, and ht quality variables are used as input to the MVA. The selection 
of variables and the training have been done separately for the different categories. For 
the MVA optimisation, the hgure of merit S/\/S + B is used, where S is the expected 
number of signal candidates in the signal region dehned by a ±1.5 MeV window around 
the known Am value 18 , while B is the number of background candidates in the signal 


region. The optimal points have signal efficiencies of 95%, 51%, 47%, and 37%, and 
background retentions of 50%, 0.33%, 0.40%, and 0.73% for LLtrig, LL, LD, and DD, 
respectively. The large background retention for the LLtrig category is due to the fact 
that the dedicated trigger intrinsically has a much lower background level. The control 
channel )■ 77“vr^ is selected by cuts on kinematic quantities, vertex quality variables, 
geometric quantities, and decay time variables. In addition, the same mass requirements 
and hducial cuts on the slow pion are applied as in case of the signal channel. Due to the 
large number of control channel candidates, which is much larger than needed for this 
analysis, 1 % of candidates are accepted at random. 


4 Asymmetry measurement 

The CP asymmetry is obtained as 


_N+ -N- 
~ N+ + N- 


( 1 ) 


for each category, where N~^ (A^~) is the yield determined from a ht to the data for a 
positive (negative) charge of the slow pion. 

The yields are determined from an extended unbinned maximum likelihood ht to the 
Am distribution. In the ht model, the signal is described by a sum of three Gaussian 
functions, where the two narrower ones are required to have the same mean value. The 
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Am [MeV] An [MeV] 

Figure 1: Distribution of Am for the control channel K~'k^ in (left) linear and in (right) 

logarithmic scale. The solid (black) line corresponds to the total fit, the dashed (blue) line 
corresponds to the background, and the dash-dotted (red) line represents the signal contribution. 


parameters of the narrowest Gaussian function and the mean value of the widest one are 
allowed to float, while the ratios between widths, as well as those between normalisations, 
of all three are constrained to the values found in the simulation. The background is 
parametrised by the product of an exponential function and a power law for the phase-space 
threshold at the pion mass 

/bg = Gbg(Am - . (2) 

Here Gbg, p and a are determined by the £t. Independent hts are performed for the four 
categories. In each category, the background parameters and the shape parameters of the 
signal component are shared between the two charges of the slow pion. 

The Am distribution for the control channel, summed over both charges of the slow 
pion with the ht function overlaid is shown in Fig.[^ Figures and show the Am 
distributions and the ht for each of the four categories and the two slow-pion charges. 
Table lists the yields from the nominal hts and the resulting asymmetries. To obtain 
the hnal result, the asymmetries of the four signal categories are combined by taking the 
weighted mean. 

5 Systematic uncertainties 

The main sources of systematic ehects are due to production and detection asymmetries 
and possible biases in the signal extraction method. The systematic uncertainty related 
to the signal extraction is estimated by comparing the nominal ht with an alternative 
one, where outside of the signal region of ±1.5 MeV around the known Am-value, only 
the background component is htted. The signal yield is obtained by subtracting the 
background extrapolated into the signal region from the total observed yield. For the 
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Figure 2: Distributions of Am split into (left) (right) D*~ and (top) LL, (bottom) LLtrig, 
including the fit function. The solid (black) line corresponds to the total fit, the dashed (blue) line 
corresponds to the background, and the dash-dotted (red) line represents the signal contribution. 


combined CP asymmetry a difference of 0.019 is found, which is assigned as a systematic 
uncertainty. 

The systematic effects that arise due to the slow pion charge asymmetry in the detector 
and a possible charge asymmetry of D* production in pp collisions in the LHCb acceptance 
are determined using the control channel. However, the control channel contains charged 
kaons which introduce an additional detection asymmetry, as the interaction cross-sections 
of and K~ with the detector material are different. In Ref. j^, the charged kaon 
detection asymmetry has been measured to be in the range 0.008 to 0.012. Assuming 
the pion detection asymmetry to be negligible, and including possible trigger effects, a 
correction of —0.010 ± 0.005 is applied to the observed asymmetry in the control channel, 
resulting in a corrected value of —0.009 ± 0.005. The absolute value of this number and 
its uncertainty are added in quadrature and assigned as a conservative estimate of the 
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Figure 3: Distributions of Am split into (left) (right) D*~ and (top) LD, (bottom) DD, 
including the fit function. The solid (black) line corresponds to the total fit, the dashed (blue) line 
corresponds to the background, while the dash-dotted (red) line represents the signal contribution. 


systematic uncertainty due to production and detection asymmetries. 

Other checks have been performed but found to have statistically insignihcant effects. 
These tests include the split into different trigger types, different run periods, and different 
magnet polarities. Also the effect of a possible difference in contamination by charm from 
beauty decays between signal and the control channel has been checked and found to 
be negligible. The total systematic uncertainty is calculated from the quadratic sum of 
the two dominant contributions, which come from the signal extraction (0.019) and the 
detection and production asymmetry (0.011), giving 0.022 for the total. 
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Table 2: Number of signal candidates and CP asymmetry obtained from the fits in the four 
categories. 


Gategory 

N+ 

N- 

•A-cp 

LL 

86± 11 

86± 12 

0.00 ±0.09 

LD 

82± 14 

83± 13 

-0.00 ±0.11 

DD 

29± 14 

66± 14 

-0.39 ±0.23 

LLtrig 

96± 11 

99± 11 

-0.02 ±0.08 

combined 



-0.029 ±0.052 


6 Result 

The time-integrated CP asymmetry in the decay )■ is determined to be 

Acp = -0.029 ± 0.052 ± 0.022 , 

where the first uncertainty is statistical and the second systematic. The result is consistent 
with no CP violation and with Standard Model expectations Q. This is the single best 
measurement of this quantity to date, with an uncertainty more than three times smaller 
than the previous determination |^. 
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